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SECTION  I 
Introduction 


The  kernel  used  in  integral  equations  can  often  represent  the  analyst's 
greatest  difficulty.  In  electromagnetics,  it  contains  a singularity  over 
which  an  integration  must  be  performed,  and  as  a result,  extreme  care  must  be 
exercised  in  performing  the  integration. 

In  this  report,  the  defining  integral  equation  in  terms  of  electric  field 
for  problems  involving  conductors  is 


n x EinC(r) 


4TTjOJ£ 


n x r 

Ja 


J(r')  • (V'V'  + k I) g (r,r ' ) 


d2r' 


(res).  (1) 


Thus,  for  a cylinder  of  radius  a with  azimuthally  symmetric  excitation  that  is 
axially  aligned  (z),  one  must  have 


^inc,  \ 
Ez  (z) 


4irjue 


r 

32  2 

(I(z') 

■'S 

2 +k 
Uz’2 

2iTa 


g (r,r ' ) d2r'  (zes) , 


(2) 


with  the  azimuthally  directed  component  of  current  density  existing  only  as 
the  solution  of  a homogeneous  equation.  In  this  expression,  k = 2tt/X,  i.e., 
the  wavenumber  in  the  ambient  medium.  Hence,  for  a circular  cylinder, 

ElnC(z)  =--i-fdz'  I(z’)  f-^+k2)  f dr  g(z,z'  <M’)  (zee),  (3) 

2 ja,£Jc  W'2  / Sir2  Jo 

with 


-jk/(z-z')2+4a2sin24>'/2 

g(z,z',  p = 0,$')  = e . =z: — • (4) 

t/(z-z  ' )2+Aa2sln2p  ' /2 

Equation  (3)  with  Eq.  (4)  is  an  exact  representation  in  this  case,  as 
there  are  no  approximations  involved  in  the  derivation,  and,  hence,  there 
are  no  restrictions  on  any  parameters  in  the  problem. 


In  this  report,  we  define  the  kernel  of  Eq.  (3)  to  be 

i r ^ 

G(z,z')  » — j | d^'gCz.z’,  <j)  = 0,4>')  d4> ' , 

8tt  J 0 

that  is,  a term  proportional  to  the  azimuthal  integral  of  the  free-space 
Green's  function.  There  is  a great  deal  of  interest  in  not  only  the 
approximations  of  the  kernel  G(z,z')  but  also  in  the  value  of  the  double 
integral  in  Eq.  (3).  However,  the  integral 


Jft  dz'  I(z’)  + k2 j G(z ,z ' ) 


has  an  analytical  solution  in  terms  of  values  at  the  end  of  the  range  of 

1 

integration  in  z'  when  I(z')  is  of  sinusoidal  (sin  kz ' or  cos  kz')  form.x 
Therefore,  we  initially  restrict  our  attention  to  the  evaluation  of  the 
integrals  in  Eq.  (3)  for  the  case  where  I(z')  is  a constant  this  gives  rise 
to  nonanalytically  integrable  integrands.  We  are  interested  in  the  value  of 
Q(z)  , i.e., 

Q(z)  =|  dz'  ^7  d4>'  g(z.z'.  * = 0,<j5’ ) . (7) 

Ja  8 it  ^0 

Because  the  major  difficulties  in  the  evaluation  of  Q(z)  occur  when  the 
singularity  in  g(z,z',  $,  <j> ' ) is  in  the  range  [a, 8],  we  focus  most  closely  on 
this  range.  However,  difficulties  have  also  been  observed  when  the 
singularity  is  outside,  but  close  to,  this  range.  We  comment  on  this 
difficulty  where  appropriate. 

Two  basic  options  exist  in  the  evaluation  of  Q(z) . Specifically,  they 
involve  the  order  i.»  which  the  $'  and  z'  integrations  in  Eq.  (7)  are  performed. 
In  this  report,  we  will  deal  with  three  different  approaches  to  the  evaluation, 
two  of  which  evaluate  the  integrals  in  the  order  indicated  in  Eq.  (7). 

Alternatives  also  exist  regarding  the  manner  in  which  g(z,  z',  j) , $')  is 
written.  One  of  the  basic  goals  of  this  report  is  the  definition  of  some 
forms  for  Q(z)  that  satisfy  certain  accuracy  requirements  over  ranges  of  the 


2 


variables  a/A,  3/A,  and  the  radius  of  the  cylinder  normalized  to  wavelength 
(a/A) . It  is  evident,  both  in  the  derivations  and  the  numerical  computations, 
that  the  particular  representations  used  for  g(z,  z',  <p,  <p')  can  affect  the 
ultimate  accuracy  in  Q(z). 

For  purposes  of  comparison,  we  derive  three  distinct  approaches  that  lead 
to  three  different  forms  of  Q(z) . Naturally,  each  approach  includes  several 
different  results,  each  of  which  is  the  outcome  of  retaining  higher-order 
terms  in  series  representations.  Also,  because  the  values  of  G(z,z')  and 
3G/3z  are  needed  in  certain  cases,  these  functions  are  considered  when 
appropriate.  For  certain  approximations,  the  results  of  numerical  computations 
are  analyzed  and  the  attendant  errors  studied. 
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SECTION  II 
Representations  for  Q(z) 


THE  CLASSICAL  APPROACH 

The  classical  approach  has  been  the  most  widely  used  and  was  one  of  the 

first  introduced  in  the  study  of  the  kernel  of  the  integral  equations  with 

which  we  are  concerned.  The  essence  of  the  method  is  described  in  Schelkunoff 
2 

and  Friis,  and  is  explicitly  carried  out  to  various  orders  of  approximation 

3 

in  Poggio  and  Mayes.  In  the  latter  work,  g(z,  z',  $,  0 ' ) as  given  by  Eq.  (4), 
is  written  in  the  form 


g(z,Zf ,0,0') 


(8) 


with 

/ 2 . 2 t , rj 

R = /4a  sm  + (z-z  ) 

Attention  is^ then  focussed  on  the  evaluation  of  the  kernel 

1 f2TT 

G(z  ,z  ' ) = — 2 I d0 ’ g ( z , z ' ,0  = 0,$').  (9) 

8-n  Jo 

when  the  last  term  in  Eq.  (8)  is  expanded  in  a Maclaurin  series  about  kR  = 0. 

A derivation  and  comparison  of  the  various  approximations  to  G(z,z')  are 
provided  in  Appendix  A and  are  given  in  Table  1. 


‘'"J.  A.  Stratton,  Electromagnetic  Theory,  (McGraw  Hill,  New  York,  1941). 

2 

S.  A.  Schelkunoff  and  H.  T.  Friis,  Antenna  Theory  and  Practice,  (John  Wiley 
and  Sons,  Inc.,  New  York,  1952). 

3 

A.  J.  Poggio,  and  P.  E.  Mayes,  Numerical  Solution  of  Integral  Equations  of 
Dipole  and  Slot  Antennas  Including  Active  and  Passive  Loading,  Air  Force 
Avionics  Laboratory  Rept.,  Wright-Patterson  Air  Force  Base,  Ohio  (1969). 
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Using  any  of  these  approximations,  one  is  still  faced  with  the  z'  integral 

indicated  in  Eq.  (7).  The  first  three  representations,  G ,,  G and  G ., 

al  az  aJ 

are  independent  of  a/A  and  the  first  term  in  the  approximation  of  G(z,z')  for 
z = z'  represents  a static  interaction.  These  three  approximations  derive 
from  the  decomposition  of  <j)(z,z',  4>  = 0,<t>'),  in  the  form  given  in  Eq.  (8), 
with  the  retention  of  only  the  first  term  of  the  MacLaurin  series  in  the 
(1  - e J )/R  term.  The  1/R  term  yields  an  exact  form,  namely,  an  elliptic 
integral,  while  the  second  term  is  0(1/ A).  The  representations  are  written  in 
an  ascending  series  in  1/An. 


AN  ALTERNATIVE  APPROACH 

The  foundation  for  the  derivation  of  an  alternative  representation  of 

4 

Q(z)  for  small  (z-z')  has  been  provided  by  Tesche.  The  approach  consists  of 
first  evaluating  the  z'  integral  in  Eq.  (7)  and  then  performing  the  <J> ' 
integration. 

In  this  representation. 


Q(z)  = ~K  ( d<D'  P(z,<t> ' ) , 
8tt  J 0 


where 


P(z  ,<p 


■>-  f 


dz'  g ( z , z ,d>  = 0,$') 


The  derivation,  detailed  in  Appendix  B,  requires  an  expansion  of  g(z,z',  <p,  $') 
about  kR  = 0 and  yields,  for  P(z,$'), 


^ r8v(<*2  sin2f  +,A 

m=0  m-  Jz-a  ' 2 ' 


The  case  where  the  singularity  is  at  the  center  of  the  range  in  E,  is  of 
particular  interest  because,  in  a collocation  solution  of  the  integral 
equation,  the  observation  points  are  usually  placed  at  the  center  of  the 

z 

F.  M.  Tesche,  Evaluation  of  the  Surface  Integral  Occurring  in  the  E-field 
Integral  Equations  for  Thin-Wire  Antennas,  Air  Force  Weapons  Laboratory, 
Kirtland  Air  Force  Base,  New  Mexico,  Mathematics  Note  29  (1973). 
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intervals.  For  this  case,  a = z - A/2  and  8 = z + A/2,  and  the  results  for 
Q(z)  are  somewhat  simplified.  We  include  below  a tabulation  of  the  representa- 
tions of  Q(z) , denoting  them  Q (z)  to  indicate  the  order  of  the  terms  retained 

in 

in  the  derivation.  Details  are  in  Appendix  B. 


V2) " ^ + s 2 + 
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,7T/2  A / 2 /A\2 
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n (Z)  = + iai  h _ its! 

V ' 4 it  2 it  I 2 


:] 


(ka)2  - (kAXka)^!  + ^ 


7 i 


,TT^2  2 2 2 A / 2 /A\2 

di^(l-k  a sin  ip)  Hn  — + a /sin  4*  + 


Q(z)=dM  + Aai  A . iV 

V J 47 r \A  3 72  / 2r  \ 2 , 


+ (ka)2  - (kA)  (ka) 

4 it2 


1 + (~f)  E/  — 

w trw 


,TT^2  2 2 2 A / 2 / A \ 2 

dip (1  - k a sin  ip)  2,n  — + Wsin  ip  + I— | 


Note  that  Q. (z)  is  independent  of  ka  but  dependent  on  kA  and  the  ratio  A/a; 

1 2 

that  Q^(z)  has  terms  of  order  ka  and  (ka)  ; and  that  the  integrand  in  Q^(z)  is 
only  slightly  more  complicated  than  that  in  Q^(z).  Furthermore,  all  the 
integrands  are  nonsingular  and  can  be  evaluated  numerically  with  ease. 
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AN  EXTENDED,  THIN-WIRE  KERNEL 

Deficiencies  have  been  observed  in  the  thin-wire  kernel,  as  illustrated 
in  Appendix  A.  There,  the  thin-wire  kernel  given  by 


-jk  /( z-z ' ) 2+a“ 

G(z  ,z  ’ ) (16) 

4tt  /(z-z  ’ )2+a2 

is  not  a suitable  representation  of  the  exact  G(z,z')  of  Eqs . (4)  and  (5)  for 
small  values  of  |z-z'|/a. 

In  our  investigation  of  representations  of  G(z,z')  or  Q(z),  we  rederived 
and  extended  the  thin-wire  kernel.  An  integral  representation  of  the  function 
g(z,z',<t>  = 0,41')  is  available  in  the  forn/ 


-jk  / (z-z ' ) 2+p2 
\4-z')2+p2 


_1_ 

25 


fdv  e^l«fl 

•/— OO 


(17) 


where 


„ <P’ 

p = 2a  sin  “ . 

The  corresponding  integral  form  for  G(z,z'),  with  an  interchange  in  the  order 
of  integration,  is  therefore 


G(z ,z ' ) 


-ij-  f dv  eJvl2-z'l 

1 6 TT  ■ j -t-oo 


H(2) 


• (18) 


The  integral  has  a well-known  result  so  that 


G(z,z ' ) 


dv  e 


jv|z-z' 


(19) 


^1.  S.  Gradshtevn  and  I.  M.  Ryzhik,  Tables  of  Integrals,  Series,  and 
Products  (Academic  Press,  New  York,  1965). 
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This  approach  to  obtain  the  exact  integral  representation  of  G(z,z')  is 
outlined  in  Hallen. 

The  substitution  of  the  series  representation  for  the  zeroeth-order  Bessel 
function,  i.e.. 


fc/k2-v2)  = ]T  (-1)1— 


2i 


i*0  22li!  F (i+1) 


2 2 

(k -v  ) , 


(20) 


allows  Eq.  (19)  to  be  written  as 


tO  i 2i  y»°0 

«.,.•) -Jj.  £ <:»  f *,.j 

8 J i=0  2 i!  m+1)  •'-*> 


i=0  2 i!  T(i+1) 


Jv  j z-z ' 1 (k2_v2ji  r(2) 
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2 (21) 


Furthermore,  one  can  introduce  the  series  representation 


„ 2 2,1  . 2i  ^ . ..2  i!  22  .-21 
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£=0 


(22) 


and  obtain 


oo 
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E (-D 


2 i! 
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-22 
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•'—00 
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(23) 


However,  the  V integral  is  nothing  more  than 


(-1) 


2 9 


22 


9z 


22 


f dv^l2-*'!  ho<2>  , 


(24) 


E.  Hallen,  Electromagnetic  Theory  (John -Wiley  and  Sons,  1962) 
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and,  using  Eq.  (17) 


G(z,z ' ) 


4tt 


£ (-D1 


i=0 


(ka)21 

22ii!  r(i+l) 


i:  ■ -2£ 

2,!  (i-E) ! 


a2* 

3z2* 


-jkXz-z ' )2+a2 
/z-z')2+a2 


This  expression  for  G(z,z')  can  be  written  in  the  compact  form 


00 

G(z,z')  = Z (-1)1 


(ka) 


2i 


i=0 


* . 

2 i!  r(i+l)  \ k / 


thin’ 


where 


’thin 


a-jki/(z-z')2+a2 

/(z-z')2+a2 


and 


3z 


2 * 


In  Eq.  (25),  the  i = 0 term  corresponds  to  the  common,  thin-wire  form 


while  retaining  the  i 
thin-wire  kernel 


V2’2')  ®thin  ’ 


0 and  i = 1 terms  corresponds  to  the  "extended 


G. (z,z*) 


gthin  ’ 


Expanding  the  terms  in  G^(z,z')  yields 


(29) 


The  extended  thin-wire  kernel  derived  above  does  not  exhibit  the 
appropriate  singularity  at  z = z*  and,  hence,  cannot  be  expected  to  yield 
suitably  accurate  results  for  Q(z)  near  that  point.  We  can  expect  however, 
that  its  range  of  validity  will  exceed  that  pertaining  to  the  thin-wire  kernel. 


THE  SMALL  DISTANCE  LIMIT  OF  THE  THIN-WIRE  FORM 

The  thin-wire  form  for  G(z,z')  was  derived  in  the  previous  section.  We 

now  consider  an  approximation  to  Q(z)  that  is  obtained  by  expanding  the 

exponential  in  a Maclaurin  series.  This  approach  is  presented  by  Harrington 

4 

and  used  for  comparison  by  Tesche. 

We  are  interested  in  the  approximation 


Q4(z) 


e~jkv/(z-z ' )^+a2 
'J  (z-z')*+a2 


(30) 


Using  the  Maclaurin  series,  one  obtains  the  integral 


<VZ> 


dz' 


V(z-z')2+a2 


(31) 


F.  Harrington,  Field  Computation  by  Moment  Methods 


(MacMillan, 


1968)  . 
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SECTION  III 

A Numerical  Comparison  of  Various  Kernels 

To  compare  the  various  kernels  or  the  functions  Q(z),  numerical  computa- 
tions are  necessary.  We  performed  these  computations  for  selected  representa- 
tions and  can  therefore  suggest  forms  that  satisfy  certain  accuracy  constraints. 
In  the  numerical  studies,  we  included  the  functions  Q(z)  shown  below. 

Table  2.  The  functions  and  their  definitions. 


Q(z) 

Eq. 

(7)  and  Eq 

. (4) (Exact) 

Vz) 

J*dz 

' G0(z,z’) 

Eq.  (27) 

Ql(z) 

Eq . 

(13) 

Q2(z) 

Eq. 

(14) 

Q3(z) 

Eq. 

(15) 

q4(*> 

Eq . 

(32) 

Qe(*> 

/ dz 

' G1(z,z') 

Eq.  (28) 

The  relative  error  is  used  as  a criterion  of  accuracy  in  the  study  of  the 
representations  and  is  defined  as 


Relative  Error 


(34) 


where  Q is  the  exact  representation  [Eqs.  (7)  and  (4)]  and  is  the  respective 
approximation.  Because  the  error  is  dependent  on  z,  a,  8,  k,  and  a,  we 
restrict  our  attention  to 

• Self  terms:  a ■ z - A/2,  B = z + A/2 

• Adjacent  terms:  a=z+A/2,  8 = z + 3A/2 

• Next-to-adjacent 

terms:  a = z + 3A/2,  8 = z + 5A/2, 


with  variations  in  wavenumber  k and  radius  a.  In  this  manner,  the  dependency 
of  Q on  z is  contained  implicitly  in  the  specific  term  we  are  considering  and 
in  the  range  [a, 8].  Therefore,  we  have  parameterized  our  curves  with  respect 
to  the  term  considered  and  to  ka,  and  have  plotted  the  relative  error  versus 
A/a. 
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RELATIVE  ERRORS  IN  THE  SELF  TERMS 

The  relative  errors  in  the  evaluation  of  the  self  terms,  Q(z)  for  a = z - 
A/2  and  8 = z + A/2,  have  been  calculated.  The  results  are  shown  in  Figs.  1 
through  8:  the  relative  errors  are  plotted  vs  A/a,  kA,  and  N/A.  The  latter 
abscissa  is  the  number  of  intervals  contained  in  one  wavelength.  Note  that 
kA  ■ (ka)(A/a)  and  N/A  = 2ir/kA.  In  Figs.  1 through  5 there  is  a region 
denoting  error  bounderies  for  certain  representations.  The  relative  error 
curve  for  each  representation  oscillated  within  these  boundaries  and  the 
detailed  nature  of  the  curves  themselves  did  not  convey  any  seemingly  useful 
information. 

The  range  of  values  of  A/a  shown  in  Figs.  1 and  2 do  not  in  reality 
correspond  to  the  values  likely  to  be  encountered  in  practice.  Because  one 
generally  uses  on  the  order  of  five  to  tens  of  intervals  per  wavelength, 
most  of  the  representations  will  suffice.  In  fact,  most  will  provide  relative 
errors  less  than  1%  in  the  commonly  used  range  of  N/A. 

For  thicker  wires,  the  relative  accuracies  begin  to  deteriorate. 

Although  the  onset  is  seen  in  Fig.  3,  it  is  more  clearly  seen  in  Figs.  4 
through  8 where  the  relative  errors  over  the  indicated  range  of  A/a  are  not 
monotonic  functions  and  often  show  a decreasing,  and  then  an  increasing, 
relative  error  as  A/a  decreases.  The  deficiency  in  is  evident  in  Figs. 

4 through  8.  Furthermore,  the  relative  inaccuracy  in  Q , Q2 , and  is 
seen  to  increase  with  increasing  ka,  and  kA  (in  general).  Inspection  of  the 
thin-wire  representation  also  indicates  a problem  as  its  relative  error 
curves  exhibit  minima.  However,  the  relative  error  in  the  extended  thin- 
wire  representation  is  reasonably  constant  with  respect  to  variations  in 
' ka  and  the  relative  error  is  less  than  1%  for  A/a  > 2. 

. i 

I A careful  study  of  Figs.  1 through  8 allows  us  to  establish  the  set  of 

. d 

representations  that  will  yield  relative  errors  less  than  1%  over  the  ranges 


10-4  < ka  < 0.4,  and  0.01  < kA  < 1.0  (Fig.  9). 
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RELATIVE  ERRORS  IN  THE  ADJACENT  AND  NEXT-TO- ADJACENT  TERMS 


The  various  representations  used  for  Q(z)  can  yield  errors  even  when 
the  interval  of  integration  does  not  contain  the  singularity.  Because  the 
representations  are  for  small  displacements,  this  may  not 

be  surprising.  However,  the  user  must  be  aware  that  difficulties  can  be 
encountered.  In  Appendix  A,  we  show  that  sizeable  errors  in  adjacent  terms 
can  arise  by  using  the  thin-wire  kernel  Gq(z,z').  We  now  investigate  a 
similar  occurrence  in  representations  of  Q(z) . In  this  portion  of  the  study 
we  include  only  those  representations  appearing  in  Fig.  9,  i.e.,  Qq,  Qe>  Q^. 

Figure  10  shows  the  range  of  relative  errors  induced  by  the  various 

approximations  as  a function  of  A/a  with  the  parameter  ka  varying  from  10  4 

_2 

to  10  . The  relative  error  is  almost  insensitive  to  the  variation. 

However,  Figs.  11  through  15  show  the  relative  errors  for  Qq  and  Q^, 
as  well  as  the  relative  constancy  of  the  relative  error  in  Qg  as  ka  increases 
The  results  presented  in  Figs.  10  through  15  allow  us  to  plot  the  range  of 
applicability  for  the  various  representations  when  evaluating  the  adjacent 
terms  (Fig.  16).  Finally,  because  we  have  considered  the  adjacent  terms, 
it  is  appropriate  to  consider  the  next-to-adjacent  terms,  a = z + 3A/2,  and 
3 = z + 5A/2,  Having  performed  a similar  study  as  that  preceding,  we  present 
here  only  .the  encapsulated  result,  namely,  the  plot  of  the  ranges  of 
applicability.  Figure  17  depicts  the  ranges  in  the  ka  - kA  space.  It  is 
evident  in  this  plot  that  the  region  is  enlarged  where  the  extended  thin-wire 
representation  provides  a relative  error  less  than  1%.  Naturally,  the  region 
of  validity  for  is  shrinking  because  it  is  a small  distance  approximation. 

The  results  of  our  investigation  of  the  various  representations  of 
Q(z)  are  summarized  as  follows: 

• If  A/a  > 10.0,  the  thin-wire  representation  can  be  used  everywhere. 

• If  A/a  > 2.0,  the  extended  thin-wire  representation  can  be  used 
everywhere . 

• If  A/a  < 2.0,  should  be  used  for  the  self  and  adjacent  terms  and 
for  the  next-to-adjacent  if  A/a  < 1.0. 

• For  1.0  <_  A/a  <_  2.0  is  used  for  the  self  and  adjacent  terms  and 

the  extended  thin-wire  for  the  next-to-adjacent  term. 


Relative  error  J2  Relative  error 


Fig.  14.  Relative  error  in  adjacent 
term  for  ka  = 0.4. 


Fig.  15.  Relative  error  in  adjacent 
term  for  ka  = 0.5. 
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REPRESENTATIONS  FOR  A SINUSOIDAL  CURRENT  DISTRIBUTION 


Several  commonly  used  basis  functions  for  conducting  wires  include 
sinusoidal  components.  Hence,  the  linear  current  density  in  Eq.  (3)  is 
written  as 


sin 

I(z  ’)  = cos  [k(z’-z.)]. 


However,  the  axially  directed  electric  field  of  such  a current  distribution 
can  be  simply  written  as^ 


Ez  ' I<2'>  f ' 


’2  31 

. “ 37'  G 


*2 


(35) 


Hence,  it  is  seen  that  one  requires  for  calculations  only  the  values  of 


G(z , z ' ) 


1 f2* 

= — 5-  / d4>  g (z , n ,<t>=0 ,0 ' ) 

n 8ir  “0 


(36) 


and 


3G(z ,z ’ ) 
3z  ’ 


1 r2lr  3 

— / dt))’  g(z,n, 0=0,0') , 


n 8tt 


(37) 


or,  in  expanded  form. 


G(z , z ' ) 


1 r/2  e-^r 

2 J d^'  r 

n 2-n  J0 


(38) 


and 


3G(z,z') 
3z ' 


--if 


n 2 it  •'o 


it/ 2 . ,,  ..  -jkr 

1ZI1  5 

r r r r 


(39) 


where 


J. .2  . 2 . 2 

r = V(z-n)  + 4a  sin  j 
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These  integrals  can  be  performed  easily  using  a numerical  scheme  as 
there  is  no  singularity.  Also,  to  circumvent  the  numerical  integration, 
one  can  use  approximate  forms.  For  instance,  for  large  |z  - z'|,r  = 

v('z-n)2  + a2. 


G(z.z') 


n 


(40) 


kA 

Fig.  18.  Representations  leading  to  kernel  evaluation  with  less  than 
1%  error. 


and 


3G(z,z ' ) 
3z' 


n 


1 jkr+1  z-n  e 
4tt  r r r 


(41) 
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For  small  |z  - z'|  the  approximations  in  Appendix  A can  be  used.  Naturally,  due 
regard  must  be  paid  in  both  cases  to  the  errors  involved  in  these  approximations 

Figure  18  shows  the  region  for  which  the  relative  error  realized  in 
approximating  G(z-z‘  = A/2)  [Eq.  (37)]  by  the  thin-wire  and  extended  thin- 
wire  representations  are  less  than  1%.  A similar  plot  is  provided  in 
Fig.  19  for  the  derivative  at  the  end  of  the  interval. 

It  can  be  concluded  that  G(z,z')  as  given  by  Eq.  (36)  can  be  evaluated 
numerically  for  A/a  < 4.0,  by  using  the  extended  thin-wire  version  for  A/a  > 
4.0,  and  the  thin-wire  version  for  A/a  > 14.0.  For  the  derivative  term, 
these  limits  become  A/a  > 4.0  and  30.0,  respectively. 


Fig.  19.  Representations  leading  to 
derivitive  evaluation  with 
less  than  1%  error. 


SECTION  IV 
Conclusions 


The  widely  used  thin-wire  representation  has  been  shown  to  be  deficient, 
particularly  when  displacements  from  the  singularity  are  small.  Alternative 
representations  have  been  provided  and  regions  of  validity  have  been 
delineated.  Constant  and  sinusoidal  representations  for  current  have  been 
considered  and  the  errors  in  evaluating  self  terms,  adjacent  terms  and 
next-to-adjacent  terms  have  been  evaluated.  Using  the  results  contained  here, 
a user  can  employ  the  form  that  is  consistent  with  his  accuracy  requirements, 
yet,  in  most  cases,  without  having  to  numerically  evaluate  the  integrals 
containing  singularities. 
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Appendix  A.  A Classical  Approach  to  the  Kernel  Approximation 


Introduction 


An  integro-dif ferential  equation  for  the  current  distribution  I(z)  on 
a cylindrical  dipole  antenna  of  length  2L  and  radius  a with  general  but 
azimuthally  symmetric  excitation  E^(z) , as  given  by  Ref.  1,  is 

I(z')  G(z,z',a)  dz ' = -jo)£  (z)  , 

where  the  kernel  G(z,a’,a),  because  of  the  assumed  symmetry,  has  the  exact 
representation 


:,z\a)  = f 
8ira  ~0 


rjka\/4  sin2  t ^ 


. . 2 V z-z 

4 sxn  + 

2 a 


The  characteristics  of  G(a,a',a)  near  its  singularity  (tp’  = 0,  z = z’)  play 
an  important  role  in  establishing  the  solution  for  I(z).  The  solution  of 
Eq.  (A-l)  or  alternative  forms  (e.g.,  Hallen's  integral  equation)  by  approx- 
imate numerical  techniques,  such  as  point  matching  (discretization  and 
8 9 

collocation) , ’ requires  (for  efficiency)  the  knowledge  of  the  valid 
approximations  to  G(z,z’,a)  that  do  not  require  aip'  integration.  Also,  the 
ranges  of  validity  of  the  approximations  must  be  known  and  observed  to  ensure 
accurate  solutions.  We  derive  here  various  representations  and  graphically 
compare  them  to  establish  their  respective  regions  of  validity.  An  error 
that  arises  when  a particular  approximation  is  used  will  be  considered. 


DERIVATION  OF  APPROXIMATIONS 


The  kernel  G(z,z',a)  can  be  written  as  the  sum  of  two  integrals:  one 

containing  a phase-stationary,  singular  integrand,  and  the  other,  a regular 
integrand : 

i I r ^ ^>1,*  r ^ i _ Q"Jkr  I 


--4  / f 

8ir  K0  r ^0 


K.  K.  Mei,  IEEE  Trans.  Ant,  and  Propag.  AP-13.  374  (1965). 

^S.  A.  Schelkunoff,  Advanced  Antenna  Theory  (John  Wiley  and  Sons,  New  York,  1952) 


where 


.2  . 2 ^ . , .2 
r =./4a  sin  + (z-z.) 

By  performing  a simple  transformation  and  substituting  a Maclaurin  series 

3 3 

about  kr  » 0 for  the  exponential  and  integrating  the  terms  to  order  k r in 
the  integrand,  the  kernel  can  be  written  as 

2 2 

G(z,z',a)  y-  2k1K(k1)  - 27Tjka  - E(k  ) _ 

8it  a 1 

2tt  3 3 , it  3 , ,.2  ^ . , 

— jk  a + -j  jk  a (z-z')  + . . . , (a-3 


where 


2a  / \ /4a2  + (z-z ' )2 


and  K(k^)  and  E(k^)  are  complete  elliptic  integrals  of  the  first  and  second 
kind,  respectively.  Equation  (A-3)  is  an  exact  expression  for  the  kernel 
when  all  the  terms  in  the  series  are  included.  However,  for  ease  of  com- 
putation , we  consider  various  approximations  for  the  kernel  for  limited 
ranges  of  the  variable  |z-z'|/a. 

An  approximation  valid  for  kr  <<  1 results  when  the  first  two  terms  of 

Eq.  (A-3)  are  retained.  This  approximation,  referred  to  as  G ,,  contains 

al 

the  result  of  the  first  integral  and  the  integral  pertaining  to  n = 1 in 

the  Maclaurin  series.  A simplification  of  G . can  be  realized  by  making 

al  10 

use  of  the  first  term  of  the  series  expansion  for  K(k^) , viz., 


K(k1)  = In  vi  - . 

2 

For  a small  |z-z'|/a,  we  have  1 - k^  s |z-z'|/2a,  so  that  a second 
approximation  is  obtained: 


-y-  2ki  *n 

8tt  a 1 


27Tjk.a|  , 


E.  Jahnke  and  F.  Emde,  Tables  of  Functions  (Dover  Publications,  New  York, 
1945). 
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and  because  k^  = 1,  a third  approximation  follows: 


Ja3 


8Ti2a 


2 2n 


8a 


Z-Z 


- 2iTjka 


(A-5) 


Let  G,  be  the  kernel  that  results  from  keeping  the  first  three  terms  of  Eq. 

D 

(A-3) , i.e.,  up  to  the  n * 2 term  in  the  Maclaurin  series.  This  approximation 
can  be  written  as 


2 

87r  a 


2k1K(k1)  - 2’rrjka  - 


,.2  2 

4k  a 


E(kx) 


(A-6) 


The  approximations  G^,  G^,  G^,  and  Gb  have  been  written  for  kr  <<  1.  From 
the  definition  of  r,  we  see  that  this  requires  ka  v4  + (z-z'/a)2  « 1.  For 
ka  « 1,  the  inequality  may  be  satisfied  for  quite  large  values  of  |z-z'|/a. 

The  kernel  can  also  be  approximated  over  the  remainder  of  the  range  of 
|z-z'|/a.  One  such  approximation,  which  is  generally  referred  to  as  the 
thin-wire  kernel,  is  obtained  by  a physical  approximation  in  which  the  source 
point  z ' is  considered  to  be  on  the  surface  of  the  cylinder  and  the  observation 
point  z to  be  on  the  axis  of  the  cylinder,  so  that  the  approximate  kernel 
can  be  written  as 


Another  approximation  can  be  derived  by  letting  |z-z'|/a  >>  4 sin  0'/2  for 
all  ip'  in  Eq.  (A-2) . Then, 


G 

e 


(A-8) 


This  kernel  is  identical  to  the  one  that  results  from  placing  both  source  and 
observation  points  on  the  axis  of  the  cylinder.  It  can  also  be  derived  from 
Eq.  (A-7)  by  requiring  that  |z-z'/a|  » 1. 
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COMPARISON  OF  APPROXIMATIONS  TO  THE  KERNEL 


2 

The  various  approximations  to  the  exact  kernel  G (normalized  to  8tt  a)  are 
plotted  as  a function  of  normalized  displacement  in  Figs.  A-l  and  A-2.  The 
ip'  integration  indicated  in  Eq.  (A-2)  was  performed  numerically  to  establish 
the  regions  in  which  the  approximations  most  accurately  represent  the  exact 

kernel.  The  figures  are  for  radius-to-wavelength  ratios  (a/X)  of  10  ^ and 

-2 

10  , respectively,  and  show  the  dependence  of  the  regions  of  validity  on 

the  dipole  thickness. 

Figures  (A-l)  and  (A-2)  show  that  a single  approximation,  adequate  for 
the  entire  range  of  |z-z'|/a,  does  not  exist.  They  also  indicate  that  G , 

clJL 

G,  , and  G _ are  excellent  approximations  for  small  values  of  the  normalized 
b aJ 

displacement  |z-z'|/a.  It  is  also  evident  that  the  accuracy  of  the  approxi- 
mations, especially  for  the  imaginary  part,  deteriorates  as  the  ratio 


Normalized  displacement  [(z-z')/a] 


Fig.  A-l.  A comparison  of  various  approximations  to  the  kernel  for  a/X  = 10 


Imaginary  part  of  G/8 


a/A  increases.  Because  it  is  necessary  to  use  different  approximations 
over  various  ranges  of  |z-z'|/a,  one  can,  in  view  of  the  simplicity  of  the 
expressions,  choose 


2 *”  H - 


z-z ' I /a  <0.3 


2 kjKCkp  - 2rrj ka 


I z-z' I /a  <3.0 


G(z,z’,a) 


-jka  \ 1+ 


| z-z' | /a  >3.0  . 


The  first  entry  is  included  because  it  has  an  analytically  integrable 
singularity  at  z = z'.  The  approximation  serves  well  for  !z-z'|/a  >6.0 
The  restriction  on  the  radius-to-wavelength  ratio  should  be  observed  when 
establishing  these  regions  of  validity. 


Exact  G 


\ 

\ G 

. \ G , G,,G 

X \ c d e 

\ |ma9'nary  sf' 


- -o.: 


-0.36 


s Imaginary  ^//  -0.38 

Ga2'Ga3'Gb'Gd'Ge"OsS^^art  . ***'''  1-0.39 


Real  part 


0.0L_i i_lJ — 

0.04  0.1 


a2  X 

Ga3 


Normalized  displacement  [(z-z")/a] 


Fig.  A-2.  A comparison  of  various  approximations  to  the  kernel  for  a/ 


Imaginary  part  of  G/8 
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A CRITICISM  OF  THE  THIN-WIRE  KERNEL 

Although  the  thin-wire  kernel  has  generally  been  used  over  the  entire 
range,  this  approximation  suffices  only  for  |z-z'|/a  <3.0.  For  normalized 
distances  less  than  3.0,  this  kernel  is  clearly  deficient  because  its  real  part 
does  not  exhibit  the  proper  singularity.  However,  because  we  are  generally 
concerned  with  its  integral  rather  than  its  functional  value,  the  thin-wire 
kernel  can  be  used  over  a wider  range  than  Figs.  A-l  and  A-2  indicate. 

It  is  interesting  to  consider  the  relative  error  that  arises  when  the 
thin-wire  kernel  is  used  in  evaluating  the  integral  on  the  right  hand  side  for 
= Zy  Figure  A-3  plots  the  relative  error  given  by 


.Z.+A/2  ~z . 

2 G , (z . , z ' ,a) dz ' - I 2 

)z-U2  d j X - 


~z  .+A/2 

I 2 G(z . , z ’ ,a)dz ' 
z .-A/ 2 2 


I J G(z . ,z  f ,a)dz  * 
Jz. -A/2  2 


vs  the  normalized  interval  half-length 
A/2a.  Also  plotted  is  the  relative 
error  for  evaluating  the  integra- 
tion over  the  adjacent  interval, 

z.  + A/2  < z.  + 3A/2.  In  the  matrix 
J ~ i 

solution  of  the  integral  equation, 
these  plots  of  the  relative  error  in 
the  diagonal  and  first-off  diagonal 
terms  in  the  coefficient  matrix. 

The  relative  error  in  the  diagonal 
and  first-off  diagonal  elements  can 
be  unacceptably  large  for  small 
intervals  (A/a<2.0).  The  is  not  a 
monotonically  decreasing  function, 
but  rather  suffers  a sharp  minimum 
when  the  interval  contains  the 
crossover  of  G^  and  G at  | z-z ' | /a=0 . 4. 
The  relative  errors  and  £^)  can  be 
reduced  bv  increasing  the  interval 


2 10.0 


A/2a 


Fig.  A-3.  Relative  errors  (tl  and 

resulting  from  the  use  of 
the  thin-wire  kernel. 
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length.  For  intervals  on  the  order  of  six  radii  (A/a^.O),  the  relative  errors 
are  less  than  27.. 

CONCLUSIONS 

Various  approximations  for  the  exact  kernel  G have  been  compared  and  a 
deficiency  in  the  thin-wire  kernel  pointed  out.  A consideration  of  the  errors 
involved  in  the  integral  of  the  kernel  leads  to  the  conclusion  that  the  thin- 
wire  kernel  should  be  used  only  for  sufficiently  large  intervals  (A/a<6.0).  A 
matrix  solution  of  the  integral  equation  for  the  dipole  using  might  not 

converge  as  the  number  of  segments  N increases  indefinitely,  but  rather  might 
begin  to  diverge  beyond  a certain  number  Nq.  To  increase  the  number  of 
segments,  other  approximations  for  the  kernel  near  the  singularity  must  be 
used . 
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Appendix  B.  An  Alternative  Approach  to  the  Kernel  Approximation 


An  alternate  representation  has  been  provided  by  Tesche.  There,  Q(z)  is 
written  as 


= _L_  f27T 

8tt2  Jq 


d<t>'  P(z, <}>'), 


(B-l) 


where 


P(z,4>') 


dz'  g ( z , z 1 , cf>  = 0,$')  , 


V2  2 a i 'y 

4a  sin  — + (z-z')“  , we  obtain 


(z,<t>')  = - X fZ  3 sin2'^-  + 52) 

H)  m-  -/z-a  1 2 ’ 


(z-3)  + 'Aa2  sin2  + (z-S)2 
- Jin  ■ — - jk(S-a) 

(z-a)  + v 4a2  sin2  + (z-a)2 


yZ-S 


- e ^ r d€R-1 

t ni  • I 

m=2  •'z-a 


The  first  two  terms  are  independent  of  ka,  as  in  the  previous  form,  and  the 
first  term  is  static.  Tesche  uses  these  two  terms  in  his  approximate 
representation,  in  which  he  sets  a = z - A/2  and  8 = z + A/2,  and  obtains 


?(*.♦')  - - 2 £n  f^+  ysin2  + 


- 2 In  sin  . ( B—  3 ) 
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Hence,  che  first-order  approximation  to  Q(z)  becomes 


Qx(z) 


^ + ^ ^ 2 + { d$’  Za  £ + 


4 it  0 


which  is  again,  as  expected,  independent  of  ka  but  dependent  on  kA. 
Let  us  now  write  Q(z)  with  higher-order  terms  included,  i.e.. 


/ 2 2 4)'  2 

. 2tt  (z-8)  + v/4a  sin  ~T  + (z-3) 

Q(z)  = -Jk  ^ - -±2  f d<0*  £n  ; 

877  0 (z-a)  + 74a2  sin2  + (z-a)2 


2tt  z-3  , 

f dc}>  ’ f dE,  R®'1  . 
0 z-a 


Thus,  including  terms  up  to  order  pertaining  to  m = 2,  we  obtain  a second 
order  approximation 


which  has  terms  of  order  ka  and  (ka)  . 

For  a = z - A/2  and  3 = z + A/2,  we  have 


The  last  term  can  be  analytically  integrated: 


dip 


where  E (tt/2 ,k * ) is  an  elliptic  integral  of  the  second  kind,  so  that 


It  is  well  known  that 


£n 


(sini^) 


dtjj 


- tt/2  ln2 


and  that 


tt/2 


/ 


sin2’^  Jin  (sini|0  d^  = ^ (1  - 2 Jin  2)  . 


Then,  we  have 


_ , , -IkA  . Zn2 

Vz)  = -Jr  + 1F 


(ka) ' 


+ ft  lka)2  - ^2  (ki)  <k*>  ft  + ft)  Zl~ 


tt/2 


+ -hr  f dU<  (1  - k2a2  sin“>)  In 


7T  + ./sin2*  + (tt)  +* 


O(kR) 


We  now  introduce  the  m = 3 term  in  Eq.  (B— 2 ) , which  can  be  written  as 


This  expression  should  be  more  than  adequate  for  most  self-term  evaluations 
The  required  integral  is  only  slightly  more  complicated  than  that  given  by 
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